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An alternate set of equations to describe the electrodynamics of superconductors at a macroscopic 
level is proposed. These equations resemble equations originally proposed by the London brothers 
but later discarded by them. Unlike the conventional London equations the alternate equations are 
relativistically covariant, and they can be understood as arising from the 'rigidity' of the superfluid 
wave function in a relativistically covariant microscopic theory. They predict that an internal 
'spontaneous' electric field exists in superconductors, and that externally applied electric fields, 
both longitudinal and transverse, are screened over a London penetration length, as magnetic fields 
are. The associated longitudinal dielectric function predicts a much steeper plasmon dispersion 
relation than the conventional theory, and a blue shift of the minimum plasmon frequency for small 
samples. It is argued that the conventional London equations lead to difficulties that are removed 
in the present theory, and that the proposed equations do not contradict any known experimental 
facts. Experimental tests are discussed. 

PACS numbers: 



I. INTRODUCTION 

It has been generally accepted that the electrodynam- 
ics of superconductors in the 'London limit' (where the 
response to electric and magnetic fields is local) is de- 
scribed by the London equationsll|,|2(. The nrs t London 
equation 



dt 



m. 



-E 



(la) 



describes the colisionless response of a conducting fluid of 
density n s , i.e. free acceleration of the superfluid carriers 
with charge e and mass m e , giving rise to the supercur- 
rent J s . The second London equation 



(lb) 



is obtained from Eq. (la) using Faraday's law and setting 
a time integration constant equal to zero, and leads to the 
Meissner effect. These equations together with Maxwell's 
equations 
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these equations are not correct, and propose an alternate 
set of equations. 

There is ample experimental evidence in favor of Eq. 
(lb), which leads to the Meissner effect. That equation 
is in fact preserved in our alternative theory. However, 
we argue that there is no experimental evidence for Eq. 
(la), even if it appears compelling on intuitive grounds. 
In fact, the London brothers themselves in their early 
work considered the possibility that Eq. (la) may not 
be valid for superconductors 0. However, because of the 
result of an experiment Q they discarded an alternative 
possibility and adopted both Eqs. (la) and Eq. (lb), 
which became known as the 'first' and 'second' London 
equations. 

It is useful to introduce the electric scalar and magnetic 
vector potentials <f> and A. The magnetic field is given 
by 



B = V x A 
and Eq. (2a) is equivalent to 

=± , ldl 
c at 



(3) 
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The magnetic vector potential is undefined to within the 
gradient of a scalar function. The gauge transformation 



A -»• A + V/ 



(5a) 



V • E = Airp 



V • B = 



(2c) 



(2d) 



are generally believed to determine the electrodynamic 
behavior of superconductors. In this paper we argue that 



10/ 
c dt 



(5b) 



leaves the electric and magnetic fields unchanged. 

The second London equation Eq. (lb) can be written 

as 



J 



-A 



(6) 
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However the right-hand-side of this equation is not gauge- 
invariant, while the left-hand side is. From the continuity 
equation, 



so that the total current J = J n + J s satisfies 



V • J 



dp 

at 



o 



(7) 



it can be seen that Eq. (6) is only valid with a choice of 
gauge that satisfies 



m e c dp 



n s e 



(8a) 



In particular, in a time-independent situation (i.e. when 
electric and magnetic fields are time-independent) the 
vector potential A in Eq. (6) is necessarily transverse, 



V • A = 



(8b) 



It is currently generally accepted that the London Eq. 
(6) is valid for superconductors in all situations, time- 
independent or not, with the gauge chosen so that Eq. 
(8b) holds (London gauge), and hence V • J = from 
Eq. (6). Together with Eq. (la) this is equivalent to 
assuming that no longitudinal electric fields (i.e. V • E ^ 
0) can exist inside superconductors, as can be seen by 
applying the divergence operator to Eq. (la). 

However, the possibility of an electrostatic field inside 
superconductors was recently suggested by the theory 
of hole superconductivity 0, IE 0]> which predicts that 
negative charge is expelled from the interior of super- 
conductors towards the surface. The conventional Lon- 
don electrodynamics is incompatible with that possibil- 
ity, however the alternate electrodynamics proposed here 
is not. We have already discussed some consequences for 
the electrostatic case in recent workj8|. 

Recently, Govaerts, Bertrand and Stenuit have dis- 
cussed an alternate Ginzburg-Landau formulation for su- 
perconductors that is relativistically covariant and has 
some common elements with the theory discussed here|9j. 
For the case of a uniform order parameter their equations 
reduce to the early London theory^ that allows for elec- 
trostatic fields within a penetration length of the surface 
of a superconductor, as our theory also does. However in 
contrast to the theory discussed here, the theory of Gor- 
vaets et al does not allow for electric charge nor electric 
fields deep in the interior of superconductors. 



II. DIFFICULTIES WITH THE LONDON 
MODEL 

Following London pj , let us assume that in addition 
to superfluid electrons there are normal electrons, giving 
rise to a normal current 



J n — (T n E 



(9) 



dJ c 2 p 
dt 4vrA| 
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with the London penetration depth Al given by 
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(10) 
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Applying the divergence operator to Eq. (10) and using 
Eqs. (2c) and (7) leads to 



(Pp 
dt 2 
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dp 
dt 



p = 



(12) 



London arguedyj that Eq. (12) leads to a rapid decay 
of any charge buildup in superconductors, given by the 
time scales 



t-l 2 = 27rcr„ ± [(27T(T n ) 2 



( c )2]1/2 



(13) 



the slower of which he estimated to be 10 sec, and 
consequently that any charge buildup inside the super- 
conductor can be ignored. He concluded from this that 
it is reasonable to assume p = inside superconductors, 
which from the continuity equation and Eq. (6) implies 
Eq. (8b). The same argument is given in Ref. [lCj . 

However, this numerical estimate is based on assuming 
for er„ a value appropriate for the normal state, and for 
Al its value near zero temperature. Instead, the tem- 
perature dependence of tr„ and Al should be considered. 
The conductivity of normal carriers a n should be propor- 
tional to the number of 'normal electrons', that goes to 
zero as the temperature approaches zero, hence a n — ► 
as T — > 0. On the other hand, the number of super- 
fluid electrons n s goes to zero as T approaches T c , hence 
Al — > oo as T — > T c . These facts lead to a strong temper- 
ature dependence of the relaxation times in Eq. (13), and 
in fact to the conclusion that long-lived charge fluctua- 
tions should exist both when T — > T c and when T — > 0. 

Eq. (12) describes a damped harmonic oscillator. The 
crossover between overdamped (at high T) and under- 
damped (low T) regimes occurs for 



2tr7„(T) 



Al(T) 



(14) 



As T — > T c , Al(T) — > co and a n (T) approaches its nor- 
mal state value at T c ; as T — > 0, Al(T) — > Al(0) and 
a n (T) — ► 0; hence condition Eq. (14) will always be sat- 
isfied at some temperature between and T c . For exam- 
ple, for a superconductor with low temperature London 
penetration depth Al ~ 200 A and normal state resistiv- 
ity p ~ lOptlcm the crossover would be at T/T c ~ 0.69; if 
Al ~ 20001 and p - 1000/xficm, at T/T c - 0.23. How- 
ever, no experimental signature of such a cross-over be- 
tween overdamped and underdamped charge oscillations 
at some temperature below T c has ever been reported for 
any superconductor. 
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For T approaching T c , the slower timescale in Eq. (13) 



is 



Tl 



47ra»(T) 



(15) 



so that for T sufficiently close to T c overdamped charge 
fluctuations should persist for arbitrarily long times. 
However, such a space charge would give rise to an elec- 
tric field in the interior of the superconductor and hence 
(due to Eq. (la)) to a current that grows arbitrarily 
large, destroying superconductivity. No such phenomena 
have ever been observed in superconductors close to T c 
to our knowledge. 

At low temperatures, Eq. (13) implies that under- 
damped charge oscillations should exist, with damping 
timescale 



1 



Tl 



2ircr n (T) 



(16) 



diverging as T — > 0. The frequency of these oscillations 



(17) 




and as T — > it approaches the plasma frequency 



c 



(18) 



Hence at zero temperature plasma oscillations in the su- 
perconducting state are predicted to exist forever, just 
as persistent currents. No such persistent charge oscilla- 
tions have been observed to our knowledge. 

Reference 0] recognized this difficulty and suggested 
that the London electrodynamic equations should only 
be assumed to be valid 'in situations in which the fields 
do not tend to build up a space charge', hence not in the 
regions close to T c and close to T — discussed above. 
On the other hand, the consequences of London equa- 
tions concerning magnetic fields and persistent currents 
are generally believed to be valid for arbitrary tempera- 
tures, and indeed experiments support this expectation. 
The fact that the implications of London's equations con- 
cerning the behavior of the charge density in supercon- 
ductors appear to have at best a limited range of validity 
is disturbing and suggests a fundamental inadequacy of 
these equations. 

A related difficulty with London's equations arises 
from consideration of the equation for the electric field. 
Taking the time derivative of Eq. (2b) and using (2a), 
(2c) and (la) yield 



1 d 2 E 
~&~dtF 



(19) 



For slowly varying electric fields and assuming no charge 
density in the superconductor 



V 2 E 
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(20) 



which implies that an electric field penetrates a distance 
Al, as a magnetic field does. Indeed, electromagnetic 
waves in superconductors penetrate a distance Al, hence 
Eq. (20) properly describes the screening of transverse 
electric fields. However Eq. (20) does not depend on the 
frequency and hence should remain valid in the static 
limit; but such a situation is incompatible with the first 
London equation (la), as it would lead to arbitrarily large 
currents. This then suggests that application of an ar- 
bitrarily small static or quasi-static electric field should 
lead to destruction of the superconducting state, which 
is not observed experimentally. 

In a normal metal these difficulties do not arise. Ap- 
plying the divergence operator to Eq. (9) and using (2c) 
and (7) yields 



dp 
dt 



-A"K(T n p 



(21) 



which predicts that any charge fluctuation in the in- 
terior of metals is screened over a very short time (~ 
10 -17 sees), and consequently that a longitudinal electric 
field cannot penetrate a normal metal. A static uniform 
electric field can penetrate a normal metal and it causes 
a finite current to flow whose magnitude is limited by the 
normal state resistivity. 

In summary, London's equations together with 
Maxwell's equations lead to unphysical predictions re- 
garding the behavior of superconductors in connection 
with space charges and electric fields, and to predictions 
that appear to be contradicted by experiment. How are 
these difficulties avoided in the conventional London pic- 
ture? By postulating that p = inside superconductors, 
and that applied static electric fields do not penetrate the 
superconductor [H I2I IToj. These are postulates that are 
completely independent of Eqs. (1) and (2), for which no 
obvious justification within London's theory exists. 



III. THE ALTERNATE EQUATIONS 

The possibility of an alternative to the conventional 
London equations is suggested by the fact that taking 
the time derivative of Eq. (6) and using Eq. (4) leads to 



dJ s n s e 2 )N 
ot m e 



(22) 



without making any additional assumptions on the gauge 
of A. Clearly, Eq. (22) is consistent with the existence of 
a static electric field in a superconductor, deriving from 
an electrostatic potential <fi, which will not generate an 
electric current, contrary to the prediction of Eq. (la). 
How it may be possible for a static electric field to exist 
in a superconductor without generating a time-dependent 
current is discussed in ref. [g- Assuming that Eq. (22) 
and Eq. (la) are equivalent, as is done in the conven- 
tional London theory, is tantamount to making an addi- 



4 



tional independent assumption, namely that no longitu- 
dinal electric field can exist inside superconductors. 

Consequently it seems natural to abandon Eq. (la) 
and explore the consequences of Eq. (22) in its full gen- 
erality. Starting from Eq. (22) , the second London equa- 
tion (lb) also follows, taking the curl and setting the 
time integration constant equal to zero as done by Lon- 
don. However to completely specify the problem we need 
further assumptions. Following the early London work0 
we take as fundamental equation Eq. (6) together with 
the condition that A obeys the Lorenz gauge: 



Js 



4nXl 



A 



V -A = -- 



dt 



(23a) 



(23b) 



These equations imply that the magnetic vector potential 
that enters into London's equation (23a) is transverse in 
a static situation as in London's case, but has a longitu- 
dinal component in a time-dependent situation. 

Application of the divergence operator to Eq. (23a), 
together with Eq. (23b) and the continuity Eq. (7) then 
leads to 



dp _ 1 dcf) 
dt ~ ~4Tr\ 2 L ~dt 



(24) 



and integration with respect to time to 

<f>(f,t) - Mr) = -4*\l(p(f,t) - Po (r)) (25) 

where <poif) and poif) are constants of integration. A 
possible choice would be 4>o = Po = 0- Instead, motivated 
by the theory of hole superconductivity p| 0, 0, @, we 
choose 



Po{r) = Po > 



(26) 



that is, a uniform positive constant in the interior of 
the superconductor. Equation (25) then implies that 
the electrostatic potential 4>{r,t) equals <j>o{r) when the 
charge density inside the superconductor is constant, uni- 
form and equal to po, hence from Maxwell's equations we 
deduce that <j>o(r) is given by 



W) = / d 
Iv 



V, 



Po 



(27) 



where the integral is over the volume of the supercon- 
ducting body. 

In summary, we propose that the macroscopic elec- 
trodynamic behavior of a superconductor is described by 
Eqs. (23) and a single positive number po, which together 
with Eq. (27) determines the integration constants in 
Eq. (25). po is a function of temperature, the particular 
material, and the dimensions and shape of the super- 
conducting body 8]. In the following we explore some 
consequences of this proposal. 



IV. ELECTROSTATICS 

For a static situation Eq. (25) is 

^r) = M?) - ^\ 2 L {p{r) - p ) 



(28) 



with <t>o{r) given by Eq. (27) [TTJ • Using Poisson's equa- 
tion we obtain for the charge density inside the super- 
conductor 



p{r) = po + AiW(r) 



(29) 



Inside and outside the superconductor the electrostatic 
potential obeys 

V 2 (0(r) - 0o (r)) = i(<KrO - Mr)) (30a) 



(r) = 



(30b) 



respectively. Furthermore we assume that no surface 
charges can exist in superconductors, hence that both 
<t> and its normal derivative dip/dn are continuous across 
the surface of the superconducting body. For given po, 
these equations have a unique solution for each value of 
the average charge density of the superconductor 



Pa 



~ J d 3 rp{r) 



(31) 



In particular, if p ave = po the solution is p(r) = po every- 
where inside the superconductor and <p(r) — 4>o(r), with 
(f>o given by Eq. (27), valid both inside and outside the 
superconductor. 

For the general case p aV e ^ Po the solution to these 
equations can be obtained numerically for any given body 
shape by the procedure discussed in ref. Q . For a spher- 
ical body an analytic solution exists, and we speculate 
that analytic solutions may exist for other shapes of high 
symmetry. Quite generally, Eq. (29) implies that for 
body dimensions much larger than the penetration depth 
the charge density is po deep in the interior of the super- 
conductor, and the potential is <fio(r). Deviations of p(r) 
from po exist within a layer of thickness Xl of the surface, 
to give rise to the given p ave . In particular, for a charge 
neutral superconductor, p ave = 0, excess negative charge 
will exist near the surface as discussed in ref. 0. 

The electrostatic field is obtained from the usual rela- 
tion 



E(f) = -V(j){r) 
and also satisfies the equation 

E(r) = E (r) + A|V 2 ^(r) 



(32) 



(33) 



with Eo{r) = —Wc/)o(r)- Deep in the interior, E(r) = 
Eo{r). Because of Eq. (22), no current is generated by 
this electrostatic field. 
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If an external electrostatic field is applied, the charge 
density will rearrange so as to screen the external field 
over a distance Xl from the surface. This is easily seen 
from the superposition principle, since the total electric 
field will be the sum of the original field and an added 
field E'(r) that satisfies 



with 



E'(r) = X 2 L V 2 E\r) 



(34) 



inside the superconductor, and approaches the value of 
the applied external field far from the superconductor. 
Equation (34) implies that the additional field E'(r) is 
screened within a distance Xl from the surface, just as 
an applied magnetic field would be screened. Quantita- 
tive results for general geometries can be obtained by the 
same procedure outlined in ref. |3| and will be discussed 
elsewhere. 

For the particular case of a spherical geometry the so- 
lution to these equations is easily obtained. The electro- 
static potential for a sphere of radius R and total charge 
q is 4>(r) — 4>{r) + 4>o(r), with 



4>(r) = 



Q sm/t(r/Ai) 



r <R 



(35a) 



_ Q sinh(R/X L ) 1 
K) ~f(R/X L ) R +LJ[ R 



r > R (35b) 



3R 



sinh(R/ Xi) 



a 



R 3 [l 



3^ 



R 2 sinh{R/X L y 



(38a) 



(38b) 



The induced charge density is easily obtained from Eq. 
(28) . Note that a dipole moment aE ext is induced on the 
sphere and that the polarizability a becomes increasingly 
reduced compared with the normal metal value a = R 3 
as the ratio of radius to penetration length decreases. For 
R >> Xl eq. (38b) yields a = (R — Xl) 3 as one would 
expect, and for R << Xl, a = R 5 /15X L . 



V. MAGNETOSTATICS 

The magnetostatics for our case is identical to the con- 
ventional case. From Eq. (23a), (2b) and (2d) it follows 
that 



B(r) = X 2 L V 2 B(f) 



(39) 



giving rise to the usual Meissner effect. The generated 
screening supercurrent is not sensitive to the presence of 
the electrostatic field. 



and 



VI. ELECTRODYNAMICS 



47T o 

Y p0 



(36a) 



(36b) 



Using Eqs. (23) and Maxwell's equations the follow- 
ing equations result for the electrodynamic behavior of 
superconductors: 



\7 2 B 



1 - I d 2 B 



(40a) 



f(x) = xcoshx — sinhx 



(36c) 



! (E-E ) = - T (E-E ) 



1 d 2 (E-E ) 
c 2 dt 2 



(40b) 



ioW = ^-^);KiJ (36d) 



V 2 J 



X 2 ' 



i d 2 j 

c 2 dt 2 



(40c) 



M r) = Q^;r>R 
r 



(36e) 



and the electric field and charge density follow from Eqs. 
(32) and (28). In the presence of a uniform applied elec- 
tric field E ext the potential is 4>(r) + <j)'(r, 6), with 

A 2 

4>'(r,0) = a^f(r/X L )E ext cos6 ; r < R (37a) 



V 2 (p - p ) 



\2 



(P - Po) + 



1 d 2 (p- Po ) 

c 2 dt 2 



(40d) 



so that all quantities obey exactly the same equation. 
Eqs. (40a, b,c) with p and E = are also obtained 
from the conventional London equations only if one im- 
poses the additional assumption that p = inside the 
superconductor. For example, the equation for the elec- 
tric field in London's case is 



b'(r,e) = (-^-r)E ext co S ; r > R 



(37b) 



y 2 E 



i 



E 



1 d 2 E 



AnVp 



(41) 
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instead of Eq. (40b). 

The simplicity of eqs. (40) derives from the fact that 
the theory is relativistically covariant Q. This is seen as 
follows. We define the current 4- vector in the usual way 

J = (J(f,t),icp(f,t)) (42) 

and the four- vector potential 

A = (A(r,t),i^r,t)) (43) 

The continuity equation sets the four-dimensional diver- 
gence of the four- vector J equal to zero, 

DivJ = (44) 

where the fourth derivative is d/d(ict), and the Lorenz 
gauge condition Eq. (23b) sets the divergence of the four- 
vector A to zero 

DivA = (45) 

Furthermore we define the four-vectors associated with 
the positive uniform charge density po and its associated 
current Jo, denoted by Jo, and the associated four- vector 
potential Aq. In the frame of reference where the su- 
perconducting body is at rest the spatial part of these 
four- vectors is zero, hence 

J = (0, icpo) (46a) 



A = (O,i0o (f)) (46b) 

in that reference frame. In any inertial reference frame, 
Aq and J are obtained by Lorentz-transforming Eq. 
(46) and are related by 

4-TT 

a 2 A = Jo (47) 

c 

with the d'Alembertian operator 

Q2 = V2 -?S < 48 > 

according to Maxwell's equations, just as the four-vectors 
J and A obey 

4-7T 

U 2 A = J. (49) 

c 

Our fundamental equation is then the relation between 
four-vectors 

n 2 (A-A ) = -^(A-A ) (50a) 

or, equivalently using Eqs. (47) and (49) 

J-Jo = --r^j(A-A ) (50b) 



which we propose to be valid in any inertial reference 
frame. In the frame of reference at rest with respect to 
the superconducting body, Jo and A have only time- like 
components, in another reference frame they will also 
have space-like components. The spatial and time-like 
parts of Eq. (50b) give rise to Eqs. (23a) and (25) re- 
spectively. 

Equations (40) can also be written in covariant form. 
Eqs. (40c) and (d) are 

2 (J-J ) = ^(J-Jo) (51a) 

and Eqs. (40a,b) 

a 2 (F-F ) = -^(F -F ) (51b) 

where F is the usual electromagnetic field tensor and 
Fo is the field tensor with entries Eq and for E and 
B respectively when expressed in the reference frame at 
rest with respect to the ions. 

VII. RELATIVISTIC COVARIANCE 

The fundamental equation (50) relates the relative 
motion of the superfluid and the positive background, 
with four-currents J and Jo, and associated vector po- 
tentials A and Ao respectively. It is a covariant rela- 
tion between four- vectors. This means it is valid in any 
inertial reference frame, with the four-vectors in differ- 
ent inertial frames related by the Lorentz transformation 
connecting the two frames. 

In contrast, the conventional London equation Eq. (6) 
with the condition Eq (8b) is not relativistically covari- 
ant, rather it is only valid in the reference frame where 
the superconducting solid is at rest. It can certainly be 
argued that the frame where the superconducting solid 
is at rest is a preferred reference frame, different from 
any other reference frame. (In fact, our covariant equa- 
tions (50) do recognize the special status of that refer- 
ence frame, as the only frame where the spatial part of 
the four-vector J is zero.) Hence it is not a priori obvi- 
ous that the electrodynamics of superconductors has to 
be describablc by relativistically covariant equations. 

However, imagine a superconductor with a flat surface 
of arbitrarily large extent, and an observer moving paral- 
lel to this surface. Since in the conventional London the- 
ory the superconductor is locally charge neutral, the state 
of motion of the observer relative to the superconductor 
is determined by its motion relative to lateral surfaces of 
the body that are arbitrarily far away. According to the 
conventional London equations this observer should be 
able to detect instantaneously any change in the position 
of these remote lateral surfaces. Since information cannot 
travel at speeds larger than the speed of light, this would 
imply that London's equations cannot be valid for times 
shorter than the time a light signal takes in traveling from 
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the observer to the lateral surface; this time increases as 
the dimensions of the body increase, and can become ar- 
bitrarily large unless one assumes there is some limiting 
value to the possible size of a superconductor. Clearly 
to describe the electrodynamics of superconductors with 
local equations whose range of validity depends on the 
dimensions of the body is not satisfactory. Our covariant 
equations avoid this difficulty. 

Note also that a relativistically covariant formulation 
does not make sense if po = 0, as assumed in the early 
London work[^. In that case, the fundamental equation 
D 2 A = — (Att/c) J makes no reference to the state of mo- 
tion of the solid and would describe the same physics 
irrespective of the relative motion of the superfluid and 
the solid, in contradiction with experiment. 



VIII. THE LONDON MOMENT 

The presence of a magnetic field in rotating 
superconductors presents another difficulty in the con- 
ventional London theory. A rotating su per conductor has 
a magnetic field in its interior given by 

2m e c _ 

B = —to (52) 



perc 

m 



with w the angular velocity. This is explained as fol- 
lows: in terms of the superfluid velocity v s the superfluid 
current is J s = n s ev s , so that Eq. (6) is 



-.4 



(53) 



Next one assumes that in the interior the superfluid ro- 
tates together with the lattice, so that at position r 

v s = Co x r (54) 

and replacing in Eq. (53) and taking the curl, Eq. (52) 
results. 

The problem with this explanation is that, as discussed 
earlier, the conventional London equations are not covari- 
ant, rather they are valid only in the rest frame of the 
superconducting body. However in writing Eq. (54) one 
is affirming the validity of London's equation in a frame 
that is not the rest frame of the solid, but is a particular 
inertial frame. To assume the validity of the equations 
with respect to one particular inertial frame that is not 
the rest frame of the solid but not in other inertial frames 
does not appear to be logically consistent and is reminis- 
cent of the old theories of the 'aether': one is stating 
that the superconductor 'drags' the 'aether' with it if it 
translates but not if it rotates. 



IX. LONDON RIGIDITY 



hole superconductivitiy|6|, the superfluid carriers are 
pairs of electrons with total spin 0. The Schrodinger 
equation for particles of spin is the non-relativistic 
limiting form of the more fundamental Klein-Gordon 
eauation|18j. It is then reasonable to expect that the 
proper microscopic theory to describe superconductivity 
should be consistent with Klein-Gordon theory. It is 
very remarkable that the London brothers in their early 
workQ, without knowledge of Cooper pairs, suggested 
the possible relevance of the Klein-Gordon equation to 
superconductivity. 

In Klein-Gordon theory, the components of the current 
four-vector J = (J(f,t),icp(f,t)) in the presence of the 
four-vector potential A = (A(f,t),i<j)(f,t)) are given in 
terms of the scalar wave function ty(r,t) bv[T^| 



J(f,t) 



P(r,t) = 



^-[**(~V--A(r,W 
2m % c 

¥(-?v"- -A(r, *))**] 
i c 



(55a) 



(55b) 



Deep in the interior of the superconductor we have (for 
superconductors of dimensions much larger than the pen- 
etration depth) 



p(f, t) = p a 



(r,t) = Mr) 



J(r, t) = A(r, t) = 



(56a) 



(56b) 



(56c) 



independent of any applied electric or magnetic fields. 
We now postulate, analogously to the conventional 
theory^, that the wave function ip(r,t) is 'rigid'. In 
terms of the four-dimensional gradient operator 



Grad = (V, - 



i d N 
cdi' 



what we mean is that the combination 



(57) 



(58) 



is unaffected by external electric and magnetic fields, as 
well as by proximity to the boundaries of the supercon- 
ductor. With = n s , the superfluid density, this 
assumption and Eqs. (55, 56) lead to 



In the conventional microscopic theory of 
superconductivity^, as well as in the theory of 



J(f,t) 



-A(f,t) 



(59a) 
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p(f, t) - po 



m e c* 



(59b) 



i.e. the four components of the fundamental equation 
Eq. (50b). 

We believe this is a compelling argument in favor of 
the form of the theory proposed here. It is true that 
for particles moving at speeds slow compared to the 
speed of light the Klein-Gordon equation reduces to the 
usual Schrodinger equation. However, by the same to- 
ken the Schrodinger equation satisfied by Cooper pairs 
can be viewed as a limiting case of the Klein-Gordon 
equation. In the conventional theory in the framework 
of non-relativistic quantum mechanics, 'rigidity' of the 
wave function leads to the second London equation. It 
would be unnatural to assume that the same argument 
cannot be extended to the superfluid wavefunction in its 
relativistic version, independent of the speed at which 
the superfluid electrons are moving. 



X. DIELECTRIC FUNCTION 

As discussed in previous sections, the electric potential 
in the interior of the superconductor satisfies 



1 d 2 



dt 2 



1 

AT 



(60a) 



while outside the superconductor the potential satisfies 
the usual wave equation 



V^- — 



i d 2 



c 2 dt 2 



= 



(60b) 



If a harmonic potential <f> ext (q,uj) is applied, the super- 
conductor responds with an induced potential <fi(q, w) re- 
lated to 4> ext by 



(61) 



and we obtain for the longitudinal dielectric function of 
the superconductor 



e s (?,w) 



u 2 + c 2 q 2 - lo 2 
c 2 q 2 — lo 2 



with 



c m e 1/2 
w p — 1 — \- A o ) 



v 47rn,e 2 ' 



(62) 



(63) 



the plasma frequency. For comparison, the dielectric 
function of the normal metal is given by the Linhardt 
dielectric functional 



£n(<7,^) 



47re 2 

1 + 



fk fk-\-q 



Zk+q — £k 



(64) 



with fk the Fermi function. 

Let us compare the behavior of the dielectric functions 
for the superconductor and the normal metal. In the 
static limit we have for the superconductor from Eq. (62) 



,(g,W-»0) = l + 



1 



c 2 q 2 



(65) 



For the normal metal, the zero frequency limit of the 
Linhardt function yields the Thomas Fermi dielectric 
function [T^ 



tTF(q) = 1 



1 



A TF q 



with 



\ 2 

TF 



4ire 2 g(e F ) 



(66a) 



(66b) 



with g(e F ) the density of states at the Fermi energy. Eqs. 
(65) and (66) imply that static external electric fields are 
screened over distances A l and Xtf for the superconduc- 
tor and the normal metal respectively. For free electrons 
we have, 



so that 



\ 2 

TF 



6mie 
e F 



3n 
2e^ 

1 3m e c 2 



2e F 



(67) 



(68) 



assuming the density of superconducting electrons n s is 
the same as that of normal electrons. Eq. (68) shows that 
the superconductor is much more 'rigid' than the nor- 
mal metal with respect to charge distortions: the energy 
cost involved in creating a charge distortion to screen 
an applied electric field is e F in the normal metal versus 
m e c 2 in the superconductor, resulting in the much longer 
screening length in the superconductor compared to the 
normal metal. 

The same rigidity is manifest in the dispersion relation 
for longitudinal charge oscillations. From the zero of the 
dielectric function Eq. (62) we obtain for the plasmon 
dispersion relation in the superconducting state 



2 2 

c q 



(69) 



Notably, this dispersion relation for longitudinal modes 
in the superconductor is identical to the one for trans- 
verse electromagnetic waves in this medium. In contrast, 
the zeros of the Linhardt dielectric function yield for the 
plasmon dispersion relation ^| 



3 o o 



(70) 



so that the plasmon dispersion relation is much steeper 
for the superconductor, since typically v F ~ 0.01c. We 
can also write Eqs. (69) and (70) as 



aJ tl 



1 



(71a) 
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obtained from (40d): 



9 

To 



\2 „2 



(71b) 



showing that low energy plasmons in superconductors re- 
quire wavelengths larger than Al according to the alter- 
nate equations, in contrast to normal metals where the 
wavelengths are of order Atf, i-e. interatomic distances. 
This again shows the enhanced 'rigidity' of the supercon- 
ductor with respect to charge fluctuations compared to 
the normal metal. 



XI. PLASMONS 

In the conventional London theory one deduces from 
Eq. (la) in the absence of normal state carriers 



dt 2 



"IP 







(72) 



upon taking the divergence on both sides and using the 
continuity equation. The solution of this equation is a 
charge oscillation with plasma frequency and arbitrary 
spatial distribution 



P(r, t) = p{f)e 



(73) 



In other words, the plasmon energy is independent of 
wavevector. This indicates that charge oscillations with 
arbitrarily short wavelength can be excited in the super- 
conductor according to London theory. Clearly this is 
unphysical, as one would not expect charge oscillations 
with wavelengths smaller than interelectronic spacings. 
Consequently one has to conclude that the 'perfect con- 
ductor' equation (la) necessarily has to break down at 
sufficiently short lengthscales. 

Experiments using EELS (electron energy loss spec- 
troscopy) have been performed on metals in the nor- 
mal state^l an d plasmon peaks have been observed, 
with plasmon dispersion relation approximately consis- 
tent with the prediction of the Linhardt dielectric func- 
tion Eq. (70). If London's theory was correct one would 
expect that in the superconducting state plasmon exci- 
tations energies should be independent of q, at least for 
values of q" 1 larger than interelectronic distances. 

However instead it is expected from BCS theory that 
plasmons below T c should be very similar to plasmons 
in the normal state[l5l lll| . This expectation is based 
on the fact that plasmon energies are several orders of 
magnitude larger than superconducting energy gaps, and 
as a consequence within BCS theory plasmons should 
be insensitive to the onset of the superconducting state. 
However no EELS experiments appear to have ever been 
performed on superconducting metals to verify this ex- 
pectation. 

In contrast, the counterpart to Eq. (72) with the al- 
ternate equations is the equation for the charge density 



o Ppl . 2 



c 2 S7 2 Ppi 



(74) 



where p p i is the difference between the charge density 
and its static value obtained from solution of Eq. (29). 
The right-hand side of this equation gives a 'rigidity' to 
charge oscillations that is absent in the London model. 
From Eq. (74) we obtain the dispersion relation Eq. (68) 
for plasmons in the superconducting state. 

Furthermore the allowed values of the wavevector q 
will be strongly constrained in small samples of dimen- 
sion comparable to the penetration depth. Consider for 
simplicity a small superconducting sphere of radius R. A 
plasma oscillation is of the form 



P P i{r,t) = p p i- 



and because of charge neutrality 



V 



d 3 rp p i{r,t) =0 



(75) 



(76) 



we obtain the condition on the wavevector 

tan(qR) = qR (77) 
The smallest wavevector satisfying this condition is 

4.493 



R 



(78) 



so that the smallest frequency plasmon has frequency 



2°4 



1 



(79) 



This shift in the plasmon frequency can be very large for 
small samples. For example, for a sphere of radius R = 
IOAl, Eq. (74) yields a 20% blue shift in the minimum 
plasmon frequency. 

The optical response of small samples will also be 
different in our theory. Electromagnetic waves excite 
surface plasmons in small metallic particles, and res- 
onance frequencies depend on sample shape and its 
polar iz ability p"?j . As the simplest example, for a spheri- 
cal sample the resonance frequency is given bv|l7| 



2 

U1 M 



Ma 



(80) 



with Q the total mobile charge, M its mass and a = i? 3 
the static polarizability of a sphere of radius R. We 
expect the polarizability to become smaller in the su- 
perconducting state as given by eq. (38b), hence our 
theory predicts an increase in surface plasmon resonance 
frequency upon entering the superconducting state for 
small samples, which should be seen for example in pho- 
toabsorption spectra. For example, for samples of radius 
IOOAl and IOAl the decrease in polarizability predicted 
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by Eq. (38b) is 6% and 27% respectively. The conven- 
tional theory would predict no such change. 

In summary, the conventional theory and our theory 
lead to very different consequences concerning the behav- 
ior of plasmon excitations when a normal metal is cooled 
into the superconducting state. In the London theory 
plasmons are predicted to be completely dispersionless. 
Within BCS theory, no change with respect to the normal 
state is expected either in the plasmon dispersion relation 
nor in the long wavelength limit of the plasmon frequency 
. Instead, in our theory the plasmon dispersion should 
be much steeper than in the normal state. Furthermore 
the minimum volume plasmon frequency should become 
larger as the sample becomes smaller, and surface plas- 
mon resonance frequencies should also become larger for 
small samples. 

XII. EXPERIMENTAL TESTS 

We do not know of any existing experiments that would 
be incompatible with the proposed theory. Here we sum- 
marize the salient features of the theory that may be 
amenable to experimental verification. 

A. Screening of applied electric field 

Our equations predict that longitudinal electric fields 
should be screened over distances Al rather than the 
much shorter Thomas Fermi length. This could be tested 
by measuring changes in capacitance of a capacitor with 
superconducting metal plates, or with a superconduc- 
tor in the region between plates, upon onset of super- 
conductivity. Such an experiment was performed by H. 
London 4] in 1936 but no change was observed. We are 
not aware of any follow-up experiment. More accurate 
experiments should be possible now. 

B. Measurement of charge inhomogeneity 

The theory predicts excess negative charge within a 
penetration depth of the surface of a superconductor and 
a deficit of negative charge in the interior. It may be 
possible to detect this charge inhomogeneity by direct 
observation, for example by electron microscopy or other 
spectroscopic tools. 

C. External electric field 

For small superconducting samples of non-spherical 
shape an electric field is predicted to exist outside the 
superconductor near the surface Q, which should be de- 
tectable by electrostatic measurements. Associated with 
it there should be a force between small superconducting 
particles leading to the formation of spherical aggregates. 



D. Internal electric field 

The predicted internal electric field is small on a mi- 
croscopic scale but extends over macroscopic distances. 
Perhaps that makes it experimentally detectable. 



E. Plasmons 

Plasmon dispersion relations should be strongly af- 
fected by the transition to superconductivity, with plas- 
mons becoming much stiffer at low temperatures. Vol- 
ume and surface plasmon frequencies should increase in 
the superconducting state for small samples. 'Small' is 
defined by the value of the ratio of a typical sample di- 
mension to Al, and effects should be detectable even for 
this ratio considerably larger than unity. EELS and op- 
tical experiments should be able to detect these changes. 



F. Polar izability 

The polarizability of small samples should be smaller 
in the superconducting than in the normal state. The 
effect should be largest at low temperatures. For sam- 
ples small compared to the penetration depth the polar- 
izability should scale as V 5 ^ 3 rather than V, with V the 
volume. 



XIII. DISCUSSION 

We have proposed a fundamental reformulation of the 
conventional London electrodynamics. The proposed 
theory is relativistically covariant and embodied in the 
single equation 

nP{A-A ) = -L(A-A ) (81) 

with A the four- vector potential, and Aq the four- vector 
potential corresponding to a uniform charge density po 
at rest in the rest frame of the superconducting body. 

Similarly as the conventional London theory pj, the 
equations proposed here can be understood as arising 
from the 'rigidity' of the microscopic wave function of 
the superfluid with respect to perturbations. However, 
in our case the relevant microscopic theory is the (rela- 
tivistically covariant) Klein-Gordon theory, appropriate 
for spin Cooper pairs, rather than the non-relativistic 
Schrodinger equation. Rigidity in this framework leads 
inescapably to the new equation (25) and hence to the 
four-dimensional Eq. (81). Furthermore, rigidity in our 
context refers to both the effect of external electric and 
magnetic fields on the superfluid wave function as well 
as to the effect on it of proximity to the surface of the 
superconducting body in the absence of external fields. 
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The constant po may be viewed as a phenomenological 
parameter arising from integration of equation (24), in- 
dependent of any microscopic theory. Instead, within the 
theory of hole superconductivity po is a positivejp&ram.- 
eter determined by the microscopic physics 0, S S 111- 
The magnitude of po does not correspond to an ionic 
positive charge but is much smaller. It originates in the 
absence of a small fraction of conduction electrons from 
the bulk which, as a consequence of the 'undressing' 7j 
associated with the transition to superconductivity, have 
moved outwards to within a London penetration depth 
of the surface. In reference |8| we estimated the excess 
negative charge near the surface for Nb to be one extra 
electron per 500, 000 atoms. For a sample of 1cm radius 
this correspond to a deficit of 1 electron per 10 11 atoms 
in the bulk, which gives rise to an electric field of order 
10 6 V/ cm near the surface. This electric field is very small 
at a microscopic level, yet it gives rise to very large po- 
tential differences between different points in the interior 
of a macroscopic sample. 

The existence of A in Eq. (81), originating in the 
positive charge po, breaks charge conjugation symme- 
try. As discussed earlier, a non-zero po is necessary for 
a meaningful relativistically covariant theory. The pre- 
diction that po is positive for all superconductors follows 
from the fundamental electron-hole asymmetry of con- 
densed matter that is the focus of the theory of hole 
superconductivity [i| 0, 0] • The fact that electron-hole 
asymmetry is a fundamental aspect of superconductivity 
is already experimentally established by the fact that the 
magnetic field of rotating superconductors always points 
in direction parallel, never antiparallel, to the mechan- 
ical angular momentum 20] . 

The electrodynamic equations proposed here describe 
only the superfluid electrons. At finite temperatures be- 
low T c there will also be a normal fluid composed of 
thermally excited quasiparticles. A two-fluid model de- 
scription of the system at finite temperatures should be 



possible and lead to interesting insights. 

The theory discussed here appears to be 'simpler' than 
the conventional London theory in that it requires fewer 
independent assumptions. It is also consistent with the 
more fundamental Klein-Gordon theory, while the con- 
ventional London theory is not, and it avoids certain diffi- 
culties of the conventional London theory. We do not be- 
lieve it contradicts any known experimental facts, except 
for the 1936 experiment by H. London|4j which to our 
knowledge has never been reproduced. Also, recent re- 
markable experiments by Tao and coworkers |2l| indicate 
that the properties of superconductors in the presence of 
strong static or quasistatic electric fields are not well un- 
derstood. The theory leads to many consequences that 
are different from the conventional theory and should be 
experimentable testable, as discussed in this paper. It 
should apply to all superconductors, with the magnitude 
of the charge-conjugation symmetry breaking parameter 
Po being largest for high temperature superconductors 0. 

Recent experiments indicate that optical properties of 
certain metals in the visible range are affected by the 
onset of superconductivity |22j. This surprising coupling 
of low and high energy physics, unexpected within con- 
ventional BCS theory, was predicted by the theory of 
hole superconductivity!^. In this paper we find that 
physical phenomena associated with longitudinal plasma 
oscillations, also a high energy phenomenon, should also 
be affected by superconductivity. Further discussion of 
the consequences of this theory and its relation with the 
microscopic physics will be given in future work. 
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